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FOREWORD 

Properties  of  the  space  charge  waves  in  a  solid  relativistic  electron  beam 
propagating  in  a  cylindrical  waveguide  is  investigated,  including  the  important 
influence  of  arbitrary  wall  impedance.  The  stability  analysis  is  carried  out 
within  the  framework  of  the  linearized  Vlasov-Maxwell  equations.  In  order  to 
examine  the  influence  of  the  axial  momentum  spread  on  the  stability  behavior,  it 
is  assumed  that  all  electrons  have  a  Lorentzian  distribution  in  the  axial 
canonical  momentum.  One  of  the  most  important  features  of  the  analysis  is  that, 
for  short  axial  wavelength  perturbations,  the  eigenfunction  can  be  described  by 
a  Bessel  function.  Moreover,  the  condition  for  zero  phase  velocity  of  the  space 
charge  wave  is  also  obtained,  in  connection  with  collective  ion  acceleration. 
Space  charge  wave  properties  for  a  dielectric  loaded  waveguide  are  also 
investigated.  For  appropriate  choice  of  dielectric  constant  e  and  thickness  of 
the  dielectric  material,  it  is  shown  that  a  strong  mode  coupling  occurs, 
exhibiting  the  growth  rate  of  instability  comparable  to  the  beam  plasma 
frequency.  The  physical  mechanism  of  instability  is  the  Cherenkov  radiation. 
This  research  was  supported  by  the  Independent  Research  Fund  at  the 
Naval  Surface  Weapons  Center. 
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I.  INTRODUCTION 

In  recent  years,  there  has  been  a  growing  interest  in  the  space 
1-4 

charge  waves  in  a  relativistic  electron  bears,  in  connection  with 

collective  ion  acceleration^  ^  and  intense  microwave  generation.2  ^ 

1-4 

For  the  most  part,  previous  theoretical  analyses  of  the  space 
charge  waves  are  limited  and  incomplete.  However,  a  complete  analysis 
of  the  space  charge  wave  is  required  to  optimize  the  collective  ion 
acceleration  and  microwave  generation.  In  this  regard,  in  the 
present  article,  we  develop  a  unified  theory  of  the  space  charge 
wave  in  a  relativistic  electron  beam,  including  the  important 
influence  of  an  arbitrary  wall  impedance  on  the  stability  behavior. 

Equilibrium  and  stability  properties  are  calculated  for  the  choice 


of  electron  distribution  function  [Eq.  (3)] 


fb(H’  pe*  V 


n^A  6(H  -  -  ymc  ) 

2it2ybm  (Pz  -  YbmBbc)2  +  A2 


where  H  is  the  energy,  P  is  the  canonical  angular  momentum,  P  is 

b  z 

2  -1/2 

the  axial  canonical  momentum,  and  y,  =  (1  -  3,  )  ,  y,  n,  ,  and  A  are 

d  b  b 

constants.  The  stability  analysis  of  the  space  charge  wave  is  carried 
out  within  the  framework  of  the  linearized  Vlasov-Maxwell  equations, 
assuming  that  v/yb  <<  1,  where  v  is  Budker's  parameter.  The  formal 
dispersion  relation  (20)  of  the  space  charge  wave  for  azimuthally 
symmetric  electromagnetic  perturbation  (9/96  =  0)  is  obtained  in  Sec.  II, 


including  the  important  influence  of  finite  wall  impedance  Z,  which  is 


generally  an  arbitrary  function  of  the  eigenf requency  ui  and  axial 
wavenumber  k.  Particularly,  we  emphasize  that  Eq.  (20)  is  derived  with 


no  a  priori  assumption  that  the  beam  radius  is  much  less  than  the  wall 
radius  Rw<  In  this  regard,  the  dispersion  relation  in  Eq.  (20)  can  be 
used  to  investigate  properties  of  the  space  charge  wave  for  a  broad 
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range  of  system  parameters. 

In  Sec.  Ill,  the  resistive  wall  instability^  with  an  arbitrary 

impedance  is  investigated  in  general,  assuming  that  all  electrons  have 

the  same  value  of  axial  canonical  momentum  (A  *=  0).  It  is  shown  that 

an  Inductive  impedance  wall  is  most  unstable.  Particularly,  the  maximum 

growth  rate  of  the  general  resistive  wall  instability  occurs  at  the 

2  2  2  2  2 

axial  wavenumber  k  satisfying  2k  Rq  =  -  Sr)*  where  the  parameter 

C  *  £r  +  iCj  is  the  root  of  =  F,  J^(x)  is  the  Vessel 

function  of  the  first  kind  of  order  i  and  F  is  the  wave  admittance  at 
the  beam  surface  [Eq.  (19)]. 

Properties  of  the  space  charge  wave  in  a  perfactly  conducting 

waveguide  is  investigated  in  Sec.  IV,  in  connection  with  intense 

microwave  generation^  ^  and  collective  ion  acceleration.^  &  It  is 

found  that  the  wave  admittance  F  in  a  perfectly  conducting  waveguide  is 

purely  capacitive,  thereby  indicating  a  stable  propagation  of  the 

space  charge  wave.  However,  for  short  axial  wavelength  perturbation 

satisfying  kR^/y^  >  10  [which  is  typical  in  the  free  electron  laser 
7  8 

application  ’  ],  the  electrostatic  eigenfunction  of  the  space  charge 

wave  can  be  accurately  represented  by  a  Bessel  function.  In  the  present 
4-6 

experiments  of  collective  ion  acceleration,  the  phase  velocity 
of  the  space  charge  wave  is  initially  required  to  be  zero,  in  order  to 
trap  and  accelerate  ions.  In  this  regard,  in  Sec.  IV,  we  also  obtain 
the  condition  for  the  phase  velocity  w/k  *=  0.  That  is  [Eq.  (42)], 


4v  /  2  i \  /  r  2  .  2  2. 

—  -  (Yb  -  l)Ur  +  k  Rq). 


We  therefore  find  from  Eq.  (42)  that  the  typical  Budker's  parameter  for 
the  zero  phase  velocity  if.  order  of  unity  in  a  practical  range  of 
physical  parameters  for  collective  ion  acceleration. 
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Space  charge  wave  properties  for  a  dielectric  loaded  waveguide 

are  investigated  in  Sec.  V,  assuming  that  the  impedance  of  the  dielectric 

material  is  purely  reactive  (i.e.,  the  dielectric  constant  is  perfectly 

real).  In  a  range  of  physical  parameters,  it  is  shown  that  the  phase 

velocity  of  the  vacuum  dielectric  waveguide  mode  is  less  than  the  beam 

9  11 

velocity,  exhibiting  possibilities  of  a  Cherenkov  radiation.  ’ 

In  fact,  a  strong  mode  coupling  between  the  vacuum  dielectric  waveguide 

mode  and  the  beam  streaming  mode  occurs  at  the  axial  wavenumber 
2  2  2  2  2 

satisfying  2k  Rq  *  ^b^i  ~  ^r^'  We  t^ere^ore  conclude  that  the  Cherenkov 

radiation  in  a  dielectric  loaded  waveguide  is  a  typical  example  of  the 

inductive  impedance  (Sec.  III).  Several  points  are  noteworthy  in  the 

analysis  of  Sec.  V.  First,  the  maximum  growth  rate  of  instability  is 

order  of  the  beam  plasma  frequency.  In  this  regard,  the  Cherenkov 

radiation  can  be  utilized  to  produce  high  power  microwave.  However, 

the  growth  rate  of  instability  decreases  substantially  with  increasing 

value  of  the  axial  momentum  spread.  Second,  the  wavelength  of  the 

microwave  radiation  generated  by  this  instability  can  be  less  than  a 

centimeter  for  a  subcentimeter  beam  radius  R^.  Finally,  we  note  that 

the  growth  rate  and  bandwidth  of  instability  increase  rapidly  as  the 

surface  of  dielectric  material  approaches  to  the  beam  surface  (RA/R  -*  1) 

u  w 
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1 1 .  THEORETICAL  MODEL 

The  equilibrium  configuration  consists  of  a  relativistic  electron 

beam  column  that  is  infinite  in  axial  extent  and  aligned  parallel  to  a 

uniform  applied  magnetic  field  The  electron  beam  radius  is 

denoted  by  R_,  and  a  finite-impedance  wall  is  located  at  radius  r  =  R  . 

0  r  w 

Cylindrical  polar  coordinates  (r,0,z)  are  introduced.  Moreover,  in  the 
present  analysis,  we  assume 


u/Yb  <<  1  »  (D 

2  2 

where  v  =  N^e  /me  is  Budker's  parameter, 

Nb  =  2 it  Jj*  dr  r  n°(r)  , 

is  the  number  of  electrons  per  unit  axial  length,  nP(r)  is  the  equilibrium 

b 

electron  density,  c  is  the  speed  of  light  in  vacuo ,  -e  and  m  are  the 
electron  charge  and  rest  mass,  respectively,  and  y^mc2  is  the 
characteristic  electron  energy  in  the  laboratory  frame.  Consistent  with 
the  low- intensity  assumption  in  Eq.  (1),  we  also  assume 


u)  ,  / o>2,  <<  1 ,  (2) 

pL  cb 

2  2 

where  w  ^  =  4-ne  h^/y^m  is  the  plasma  frequency-squared  and  = 
eB0^bmC  is  electron  cyclotron  frequency. 

In  the  present  analysis,  we  investigate  stability  properties  for 
the  choice  of  equilibrium  distribution  function 


n  n  A  6(11  -  -  ymc  ) 

f°(H,Pft,P  J  =  - — -  ,  (3) 

2tt  V>  (pz  "  Vfcbc)  +  A 

2  A  2  2  T  /  2 

where  H=(mc  +  c  js  )  -  ei^(r)  is  the  total  energy,  =  r[p^  -  (e/2c)rBQ] 

is  the  canonical  angular  momentum,  pz  =  Pz  ~  (e/c)A®(r)  is  the  axial 

<5 

canonical  momentum,  $^(r)  and  Az(r)  are  the  equilibrium  electrostatic 
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and  axial  component  of  vector  potentials,  respectively,  y  =  (1  -  g^-1/2 

b  b 

Y»  n^,  and  A  are  constants. 

In  the  subsequent  perturbation  analysis,  use  is  made  of  the  linearized 
Vlasov-Maxwell  equations  for  azimuthally  symmetric  perturbations  (3/36  *  0) 
about  a  solid  electron  beam  described  by  Eq.  (3).  We  adopt  a  normal¬ 
mode  approach  in  which  all  perturbations  are  assumed  to  vary  according  to 

6\p(x,t)  =  <Kr)exp[i(kz  -  mt)]  , 


where  Imui  >  0.  Here,  ui  is  the  complex  eigenf requency  and  k  is  the  axial 
wavenumber  .  The  Maxwell  equations  for  the  perturbed  electric  and 
magnetic  field  amplitudes  can  be  expressed  as 


l  x  W  =  i(<o/c)£(jc)  , 

V  x  (l/u)£(;£)  *  (Wc)^)  -  i(w/c)c|(^c)  , 


(4) 


where  e  and  p  are  the  dielectric  constant  and  permeability,  respectively, 
|(^c)  and  J^(^)  are  the  perturbed  electric  and  magnetic  fields,  and 


W 


=  -e 


fd'3p  x 


(5) 


is  the  perturbed  current  density.  Note  that  c  =  p  =  1  in  vacuo. 
In  Eq.  (5), 


f° 

v 1  x  b(x5)1 

Vs-tf  ‘ e 

dr  exp(-iwr) 

E(x  )  +  - 

—  oo 

. 

is  the  perturbed  distribution  function,  t  *  t’  -  t,  and  the  particle 
trajectories  ^c'(t')  and  j^'(t')  satisfy  dx'/dt'  =  v'  and  dj^'/dt'  «  -ev'  * 


with  "initial"  conditions  x' (t '  =  t)  =  x  and  v’(t'  =  t)  =  v. 

'V  %  %  'V, 
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In  general,  the  permeability  p  of  a  dielectric  material  even  in  the 
wall  differs  from  unity  by  only  a  few  parts  in  10"*.  Therefore,  we 
approximate  p  ■  1  in  the  remainder  of  this  paper.  Making  use  of 
Eq.  (A),  it  is  straightforward  to  show  that 


E  (r)  =  (kc/we ) B  (r )  , 
r  u 


(7  ) 


B  (r)  *=  i[u>e/c(<A/c2  -  k2)][0E  (r)/3r]  , 
o  z 


and 


(r  h  r  It  +  *7  c  '  r2)  Ez(r)  =  4l,lk 


P (r )  -  Jt(r)j 

c  k  ' 


(8  ) 


where  is  the  azimuthal  component  of  the  perturbed  magnetic  field, 
and  and  E^  are  the  radial  and  axial  components,  respectively, 
of  the  perturbed  electric  field,  p(r)  is  the  perturbed  charge  density 
and  ^(r)  is  the  axial  component  of  the  perturbed  current  density. 

To  lowest  order,  the  axial  motion  of  the  particle  orbit  is  free- 
streaming 


z' 


(f  -  t)  . 


(9  ) 


Moreover,  within  the  context  of  Eq.  (2),  we  neglect  the  terms  proportional 
to  j£j,(r)  on  the  right-hand  side  of  Eq.  (6),  where  Ex(r)  is  the  transverse 
component  of  the  perturbed  electric  field.  Finally  we  assume  a  slow 
rotational  equilibrium  characterized  by 

(10) 


thereby 


approximating 


12 


(11) 
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in  the  arguments  of  the  perturbation  amplitudes  on  the  right-hand  side 
of  Eq.  <6). 

Substituting  Eqs.  (9)  and  (n)  into  Eq.  (s) ,  we  obtain  the  perturbed 
distribution  function 


ei  -L  ,  .  b 

f.  (r,p)  ■  - r - 7 —  E  (r)  - -  . 

b  K  w  -  kp  /vm  z  3  p 
z  rz 

Carrying  out  the  momentum  integration  with  Eqs.  (3)  and  (12),  the 
differential  equation  (8)  can  be  expressed  as 


(12) 


(r  h r  It  +  1,2  = 0  *  0  -< r 


<  k  . 


(r  h r  k+  p2K(r)  =  0  ’  R0  '  r  '  R„  * 


where 


H) 


-  kBfac  +  ikA/y^m)2 


-  1 


E  (R  ) 


Z(w,k)  =  -i(uiR  /c)  — 


z  w 


Be(Rw) 


(13) 


(14) 


2  w  .2 
P  -  —  -  k  , 

c 

and  the  beam  radius  Rq  is  defined  by 

R0  =  2c2(^  '  V'V'V’cb  “  %  -  ^pb/2yb)  •  (l5) 

For  the  convenience  in  the  future  analysis,  we  define  the 
wave  impedance  Z(u,k)  of  the  wall  as 


(16) 


where  R  is  the  radius  of  the  finite  impedance  wall.  Evidently,  the 


solutions  to  Lq.  ( 1 3 )  are  given  by 


Ez(r) 


A J  0 ( Tr )  ’  0  '  r  <  RQ  , 


BJ0(pr)  +  CNQ(pr)  ,  Rq  -  r  _<  Ry  , 


(17) 
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where  J^(x)  and  N^(x)  are  Bessel  functions  of  the  first  and  second  kind, 
respectively,  of  order  l.  Making  use  of  Eqs.  (7)  and  (16),  and  the 
boundary  conditions  of  E^(r)  at  r  *  R^,  we  obtain  the  dispersion  relation, 

tJ,(0 

1  T~(0  *  F^w,k^  »  (is) 


where  F(w,k)  is  the  wave  admittance  at  the  beam  surface  defined  by 

+  gCpR^CpR,,) 


F(uj,U)  =  pR 


0  J0(PR0)  +  g(pRw)H0(pR0)  * 


(  19) 


and 


£  = 


TR, 


(20) 


In  Eq.  (19), 


g(x) 


Jx(x) 

xNQ(x)  -  ZNx(x) 


xJQ(x) 
J,(x)  ' 

±  J 


(21) 


and  p  is  defined  in  Eq.  (14).  In  the  remainder  of  this  article,  we 
make  use  of  Eq.  (18)  to  investigate  properties  of  the  space  charge 
wave  for  various  values  of  wall  impedance  Z  in  Eq.  (16). 
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III.  INFLUENCE  OF  WAVE  ADMITTANCE  ON  THE  SPACE  CHARGK  WAVES 

For  convenience  in  the  subsequent  analysis,  we  denote  the  root  of 
Eq.  (18)  by 


5  -  Cr(oa,k)  +  U^ui.k)  , 


(22) 


for  a  specified  value  of  wave  admittance 


F  «  fexp (i$)  =  Ff  +  iFi  , 


(23  ) 


where  £r  and  £  are  the  real  and  imaginary  parts  of  £,  and  f  and  <J>  are 
the  magnitude  and  phase  angle,  respectively,  of  the  general  wave 
admittance  F  at  the  beam  surface.  From  Eq.  (18)  with  Eq.  (23),  we  note 


*  * 

£  JX(C  ) 

J0(C*) 


*  f  exp  (— i 4= )  » 


(24) 


*  * 

where  £  is  the  complex  conjugate  of  £  defined  by  £  =  -  iC^. 

Shown  in  Fig.  1  are  contours  of  constant  phase  angle  4>  and  modulus 
f  of  the  function  F  «  £J^(£)/Jq(£)  in  the  complex  plane  £  *  £f  +  i£^. 

We  note  from  Fig.  1  that  the  root  £  of  Eq.  (18)  approaches  zero  or 
as  the  magnitude  of  the  wave  admittance  is  reduced  to  zero,  f  -+  0. 

On  the  other  hand,  £  approaches  8Qn  as  f  increases  to  infinity. 

Here  ciQn  and  are  the  nth  root  of  J^u^)  =  0  and  = 

respectively.  For  specified  values  of  f  and  we  note  that  the  root 
£  -  £f  +  i£^  can  be  determined  from  Fig.  1. 

Within  the  context  of  Eq.  (1),  it  is  very  useful  in  the  subsequent 
analysis  to  note  that 


I gj  -  |  <<  kc  . 


(25) 
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In  the  remainder  of  this  section,  we  investigate  the  resistive  wall 
instability^  in  general,  assuming  that  all  electrons  have  the  same  value 
of  axial  camonical  momentum,  i.e., 


A  =  0  . 


(26) 


Making  use  of  Eqs.  (14),  ( 20) ,  and  (26).  obtain. 


n2  = 
1 


2[(y 


+  6)  +  1]  l 


y|((y  +  A)2  +  l]1//2  -  y (y  +  i)j  , 

( 27) 


for  the  resistive  wall  instability  in  general.  In  Eq.  (27),  Jl  is  the 
normalized  growth  rate  defined  by 


lm(ui  -  kB,  c) 

a.  - 7 — - — 

1  >/yb 

and  the  parameters  y  and  6  are  given  by 


(28) 


and 


y  -  , 


«  -  U2  -  C2)/2^rCi  , 


(29) 


(30) 


respectively.  In  obtaining  Eq.  (27),  use  has  been  made  of  Eq.  (25). 

2  2 

(a)  Inductive  impedance  with  Shown  in  Fig.  2  are  plots  of 

the  normalized  growth  rate  versus  y  obtained  from  Eq.  (27)  for 
2  2 

>  £r  and  several  different  values  of  6.  As  shown  in  Fig.  2, 

is  the  maximum  growth  rate  for  a  given  value  of  parameter  6,  and  y* 

m 

is  the  corresponding  value  of  y.  Obviously  from  Fig.  2,  the  maximum 

growth  rate  increases  rapidly  with  the  increasing  value  of  { 6 |  in  the 

inductive  impedance  case.  In  Fig.  3,  we  plot  the  parameter  y^ 

m 

corresponding  to  the  maximum  growth  rate  versus  6  (solid  line) 
obtained  from  Eq.  (27)  and  y^  =  -6  (dashed  line).  Evidently,  we  note 
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from  Fig.  3  that  the  maximum  growth  rate  ft™  and  the  corresponding  y 
can  be  approximated  by 


ft®  =  (|6|/2)1/2  ,  (31) 

and 

*  -*  •  (32) 


for  | 6 |  >1.  Substituting  Eqs.  (29)  and  (30)  into  Eq.  (32),  we 

conclude  for  the  inductive  impedance  that  the  maximum  growth  rate  of 

the  resistive  wall  instability  occurs  at  the  axial  wavenumber  k 

satisfying  k2R2  “  Yb^i  ~  ^r^  ^or  ^  *  ** 

2  2 

(b)  Capacitive  impedance  with  £r  >  Presented  in  Fig.  4  are 

plots  of  the  normalized  growth  rate  versus  y  obtained  from  Eq.  (27) 

2  2 

for  £  0^  and  several  values  of  parameter  6.  As  shown  in  Fig.  4, 

we  again  define  the  maximum  growth  rate  ft™  and  corresponding  y^  for  a 
given  value  of  j.  In  order  to  illustrate  the  stability  dependence 

on  the  parameter  6,  in  Fig.  5,  we  plot  the  maximum  growth  rate  ft™ 

c  2  2 

and  corresponding  y  versus  6  obtained  from  Eq.  (27)  for  £  >  The 

m  r  l 

growth  rate  increases  rapidly  with  the  decreasing  value  of  1 6  j  for  the 
capacitive  impedance  case. 

2  2 

(c)  Perfectly  resistive  with  {,  =  For  the  perfactly  resistive 

~  ~  / 

impedance  case  characterized  by 


2 

r 


(33) 


Eq.  (27)  can  be  simplified  as 


Qi - 21 — Uyl  <y2  +  D1/2  -  y2J  .  (34) 

1  2 (y^  +  1) 

which  gives  the  maximum  growth  rate 


ft®  =  0.3530  ,  (35  ) 
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at  y  ■  0.577.  We  therefore  conclude  that  the  induct've  impedance 
case  is  most  unstable.  In  the  following  sections,  we  investigate 
properties  of  the  space  charge  waves  in  perfectly  capacitive  and 
Inductive  impedances. 
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IV .  SPACE  CHARGE  WAVES  IN  A  PERFECTLY  CONDUCTING  WAVEGUIDE 

In  this  section,  we  investigate  pronerties  of  the  space  charge 
waves  in  a  relativistic  electron  beam  assuming  that  the  wall  impedance 
is  zero,  i. e. , 


Z  =  u  ,  (36) 

and  that  all  electrons  have  the  same  value  of  axial  canonical  momentum 
(A  *=  0).  Substituting  Eqs.  (21)  and  (36)  into  Eq.  (19)  and  defining 

2  „  2  2,2.2  ,  . 

q  =(k  -  (d  /c  )Rq  ,  (37 ) 

we  obtain 

I1(q)K0(Rwq/R0)  +  ^(R^/Rq)^  (q) 

F(q)  =  q  *  (38) 

where  I^(x)  and  K^(x)  are  the  modified  Bessel  functions  of  the  first 
and  second  kind,  respectively,  of  order  8..  For  a  real  value  of  q,  the 
wave  admittance  F(q)  in  Eq.  (38)  is  a  positive  real  value  (F^  =  0), 
thereby  giving  the  phase  angle  ip  *  0  in  Eq.  (23)-  In  this  regard,  the 
space  charge  wave  admittance  F  in  a  perfectly  conducting  waveguide  is  a 
perfactly  capacitive,  corresponding  to  the  horizontal  line  ^  *  0  in 
Fig.  1.  Without  further  analysis,  we  therefore  conclude  from  Fig.  4 
that  the  space  charge  mode  in  a  perfectly  conducting  waveguide  is 
stable  [6  *  <*>  in  Eq.  (  27)  J . 

Figure  6  shows  plots  of  the  admittance  F  =  Fr  versus  4 
(dashed  curves)  obtained  from  Eq.  (18)  and  Fr  versus  q  (solid  curves) 

obtained  from  Eq.  (38)  for  Rn/R  *  0.1,  0.5  and  0.9.  The  horizontal 

U  w 

scale  in  Fig.  6  represents  both  and  q.  In  Fig.  6,  n  =  1,  2,  and  3 
denote  tne  radial  mode  number  of  the  space  charge  wave.  Note  that 
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for  specified  values  of  q  and  R^/Rp,  the  parameter  £  is  determined  from 

Fig.  6.  Shown  in  Fig.  7  is  plot  of  versus  q  determined  from  Fig.  6 

f°r  Rq/rw  *  Obviously,  the  root  £r  approaches  6qr  as  the  parameter 

q  increases  to  infinity.  This  is  similar  to  Fig.  1,  since  the  wave 

admittance  F^  is  monotonically  increasing  with  the  increasing  value  of  q. 

In  order  to  illustrate  the  influence  of  the  parameter  q  on  the 

space  charge  mode,  we  present  in  Fig.  8  plots  of  the  eigenfunction  Ez(r) 

versus  r/RQ  obtained  from  Eq.  (17)  for  (a)  n  •  1,  (b>  n  *  2,  R^/R^  =0.5 

and  several  values  of  q.  The  perturbed  axial  electric  field  E  (r)  in 

z 

Fig.  8  is  normalized  by 

,  fR 

~2  Wdr  r  E  (r)  =  1. 

Ro  J° 

It  is  obvious  from  Fig.  8  that  the  axial  electric  field  E  (r)  in  Eq.  (17) 

z 

can  be  approximated  by 

(aj0(£s  r/R. )  •  0  <  r  <  R  , 

Ez(r)  -l  °  °n  0  0  (39) 

10  ,  otherwise  , 


for  q  >  10.  We  therefore  conclude  that  inside  the  beam,  the  electro¬ 
static  eigenfunction  of  the  large  axial  wavenumber  perturbations  is 
represented  by  a  Bessel  function.  For  example,  Eq.  (39)  is  an  excellent 

approximation  for  the  electrostatic  eigenfunction  of  the  free  electron 
7  8 

laser  instability.  ’ 

Making  use  of  Eq.  (20)  and  (37),  Eq.  (14)  can  be  expressed  as 

2  2 


7  W  iQ 

u-  k^c)  -  ~2  22~  —  r  • 

YbUr  q  } 

where  the  root  £r(q)  is  determined  from  Fig.  7  for  given  values  of 
Rq/R^  and  n,  and  use  has  been  made  of  L  =  0.  The  eigenf requency  u 
and  axial  wavenumber  k  are  obtained  from  the  simultaneous  solution 


(40) 


of 
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Eqs.  (37)  and  (40).  Shown  in  Fig.  9  is  a  plot  of  the  dispersion  curve 
in  the  (w,k)  parameter  space  for  n  ■  1,  R./R  ■=  0.5,  y,  »  3,  and 

U  W  D 

2  2  2 

v/y  -  id  R  /4c  «  0.1.  Simultaneously  solving  Eqs.  (37)  and  (40)  for 

D  pb  U 

q  and  £^(q)  give  two  distinct  dispersion  curves.  The  fast  wave 

mode  in  Fig.  9  corresponds  to  the  phase  velocity  V  ^  *  w/k  >  B^c 

and  the  slow  wave  mode  corresponds  to  V  ,  <  B,  c. 

ph  b 

After  some  straightforward  algebraic  manipulation  with  the 
2  2  2 

definition  v/y,  =  w  ,  R„/4c  we  can  rewrite  Eq.  (40)  by 
b  pb  u 


/iii —  e  \ 

\  kc  )  3 


i  2  ..  2  2 

1  -  u)  /k  c 


4v  _ _ 

.2  ,  . .  2,, 2  2,  ,2  2  ’ 

yb  +  (1  -  W  /k  c  )k  R0 


(41 ) 


which  determines  the  phase  velocity  V  ,  =  oi/k  of  the  space  charge 

ph 

wave  in  terms  of  the  normalized  axial  wavenumber  kR^,  the  Budker's 
parameter  v  and  energy  yb  of  the  electron  beam.  In  recent  years, 
there  has  been  a  considerable  increase  in  interest  in  collective  ion 
acceleration  by  a  slow  space-charge  wave-3  in  a  relativistic  electron 
beam.  However,  the  phase  velocity  of  the  space  charge  wave  is  initially 
required  to  be  zero,  in  order  to  trap  and  accelerate  ions.  In  this  regard, 

we  derive  the  condition  for  V  .  «  u/k  =  0  from  En.  (41).  That  is, 

.  Ph 


rh-  ^-l)^r(kR0)  +  k2Ro]  • 


(42) 


where  the  root  £  is  determined  in  terms  of  the  normalized  axial  wave- 
r 

number  kR„.  For  given  values  of  R„/R  and  n,  the  root  (,  required  in 
0  U  w  r 

Eq.  (42)  can  be  found  from  Fig.  7  where  the  horizontal  scale  q  is 
replaced  by  kR^. 

Figure  10  is  plots  of  £  versus  R  /R  for  n  *  1  and  several  values 
°  X  0  W 


of  kR^.  The  root  £r  is  monotonical ly  increasing  to  Cqj 


2.4  as  the 


ratio  R_/R  increases  from  tero  to  unity.  Making  use  of  Fig.  10, 
0  w 


25 


S'-  ... 


NSWC  TR  81-65 


the  necessary  value  of  Budker's  parameter  for  the  zero  phase  velocity 
of  space  charge  wave  is  found  from  Eq.  (42).  For  example,  for  kR^  *  1, 
Rq/Rw  m  0.5  and  y^  *  1.2,  we  find  Av/y^  «=  1.46  for  V  ^  *  w/k  ■  0. 

Typical  Budker's  parameter  for  the  zero  phase  velocity  is  order  of  unity 
in  a  practical  range  of  physical  parameters  for  collective  ion  acceleration. 
Finally,  for  a  small  phase  velocity  satisfying  |u/kc|  <<  6^ ,  Eq.  (41) 
can  be  approximated  by 


4v  .2  .  v 

—  =  (y.  -  1) 

Yk  b 


(43) 


which  is  identical  to  Eq.  (42).  However,  from  Eq .  (43),  the  Budker's 
parameter  is  determined  in  terms  of  the  oscillation  frequency  oj  and  the 


phase  velocity  V  ^  »  <u/k  required  to  ion  accelerations. 
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V.  SPACE  CHARGE  WAVES  FOR  A  DIELECTRIC  LOADED  WAVEGUIDE 


In  this  section,  we  investigate  properties  of  the  space  charge 

waves  for  a  solid  electron  beam  in  a  dielectric  loaded  waveguide. 

The  solid  electron  beam  with  radius  Rq  propagates  through  a  cylindrical 

waveguide  loaded  with  dielectric  material  in  the  range  R  <  r  <  R  . 

w  c 

A  grounded  cylindrical  conducting  wall  is  located  at  radius  R  . 

c 

Vie  approximate  the  permeability  of  a  dielectric  material  p  »  1.  In 
this  regard,  the  perturbed  axial  electric  field  Ez(r)  and  azimuthal 

magnetic  field  B  (r)  are  continuous  across  the  dielectric  boundary  at 

6 

r  -  R^.  From  Eq.  (8),  we  obtain 


(i 


3_ 
3  r 


d  r 


(r)  =  0 


(44) 


inside  the  dielectric  material  (R^  '  r  r  Rc).  l!ere  P2  is  defined  by 


2  2.2  .2 
p  =  0.  t/c  -  k 


(45) 


The  solution  to  Eq .  (44)  can  be  expressed  as 


Ez(r)  =  A[JQ(p2r) 


JQ(p2Rc)Nl(p2r)/K0(P2Rc)]  * 


(46  ) 


where  A  is  an  arbitrary  constant.  Substituting  Eq .  (46)  into  Eq.  (7), 
and  making  use  of  the  boundary  conditions  at  r  =  R^  and  Eq .  (16),  we 
obtain  the  impedance 


n  -  v^vy vn) 

t  J1(n)N0(nRc/Rw)  -  JQ(nRc/Rw)N1(n)  * 


(47  ) 


at  the  surface  of  dielectric  material  (r  =  R  ) .  In  Eq.  (47),  the 

w 

parameter  n  is  defined  by 


2 

n 


2  2 
P0R 
2  w 


2  2 
(cA/c 


k2) 


(48  ) 
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Equation  (18),  when  combined  with  Eqs.  (19),  (21),  and  (47),  yields  a 
closed  dispersion  relation  for  the  space  charge  wave  mode  in  a  dielectric 
loaded  waveguide. 

Making  use  of  Eqs.  (14),  (19),  and  (20),  it  can  be  shown  that 

2 

in  the  vacuum  dielectric  waveguide  characterized  by  in  '  0  in  Eq.  (14), 

pb 

the  dispersion  relation  in  Eq .  (18)  is  simplified  as 

g(pR  )  -  0  ,  (49) 

W 

22  2  1/2 

where  p  “  (ui  /c  -  k  )  .  Substituting  Eqs.  (21)  and  (47)  into  Eq.  (49), 

we  obtain  the  vacuum  transverse  magnetic  (TM)  mode  dispersion  relation, 

pRwJo(pV  „  V'>>V"r'c/V  -  VVW") 

— c  -l  — . . —  .  .  —  ■  - - - - -  (  ^0) 

J1  (PRw)  £  "  J0("Rc/Rw)Nl(,|)  ’ 

where  the  parameter  pR^  is  defined  by 

2  2  2  2  2  2 

PV  -  R‘(<//c  -  kZ)  .  (51) 

w  u 


It  is  instructive  to  examine  Eq .  (50)  in  the  limit  t  -*  1 . 

Making  use  of  PR  =  q»  obtain 
w 

JQ(nRc/Rw)  =  0  ,  (52) 


from  Eq .  (50)  for  e  ■  1.  Equation  (52)  gives  the  familiar  vacuum  TM 

mode  dispersion  relation. 


c 


in  a  perfectly  conducting  waveguide. 

the  limit  of  both  R  -*  R  and  R  -+  0, 
w  c  v; 

relation  in  Eq.  (53)  and, 


(53) 


Moreover,  we  can  show  that  in 
Eq .  (50)  gives  the  dispersion 


2 

pi 

2 


c 


t 


(54) 
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respectively.  Note  that  the  case  R  -►  0  corresponds  to  a  completely 
filled  dielectric  waveguide. 

For  given  values  of  the  dielectric  constant  t  and  the  ratio  R  /R  , 

w  c 

the  parameter  pR^  is  determined  from  Eq.  (50)  in  terms  of  n •  The 

oscillation  frequency  u  and  axial  wavenumber  k  in  a  vacuum  dielectric 

loaded  waveguide  are  obtained  from  the  simultaneous  solution  of 

Eqs.  (48)  and  (51)  for  specified  n  and  pR  .  Figure  11  is  plots  of  the 

w 

vacuum  HI  mode  dispersion  relation  in  the  (u,k)  parameter  space  for 

(a)  R  /R  -0.8  and  several  values  of  the  dielectric  constant  c ,  and 

w  c 

(b)  e  “  4  and  several  values  of  the  ratio  R  /R  .  Vie  remind  the  reader 

w  c 

that  the  thickness  of  dielectric  material  increases  from  aero  to  Rc 
as  the  ratio  Rw/R(,  decreases  from  unity  to  zero.  All  plots  in  Fig.  11 
correspond  to  the  first  radial  mode  number  (  n  =  1).  The  dispersion 
curve  for  R^/R^  *  0  in  Fig.  11(b)  represents  the  dispersion  relation 
of  a  completely  filled  dielectric  waveguide  [Eq.  (54)].  On  the  other 
hand,  the  curves  for  c  *  1  in  Fig.  11(a)  and  for  P  /R  =  1  in  Fig.  11(b) 

W  Q 

correspond  to  the  ordinary  dispersion  relation  in  Eq .  (53)  where  the 

phase  velocity  V  ^  =  u/k  is  always  faster  than  the  speed  of  light  (u/k  > 

However,  the  phase  velocity  of  the  dispersion  relation  in  a  dielectric 

loaded  waveguide  is  sometimes  less  than  the  speed  of  light  (u/k  <  c). 

For  example,  for  R^/R^  =0.8  and  e  =  4  in  Fig.  11(a),  we  find  u/k  >  c 

for  kRc  <3.3  and  u/k  <  c  for  kR  >  3.3. 

c 

It  is  evident  from  Eq.  (14)  and  (20)  that  in  the  vacuum  dielectric 

2 

waveguide  (u^  *  0),  we  can  equivalently  express  the  combination 
of  Eqs.  (48),  (50 ),  and  (51)  as, 

2 


2  -  k2  -  J[2(u,k)/R2  , 


(55) 
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where  £(w,k)  is  determined  from  Eq.  (18)  in  terms  of  the  wave 

admittance  F  at  r  ■  Rq.  Equation  (55)  is  a  compact  form  of 

the  vacuum  HI  mode  dispersion  relation  in  a  dielectric  loaded 

22  222  1/2 

waveguide.  Shown  in  Fig.  12  is  a  plot  of  u  »  (k  c  +  c  /R  ) 

Un  c 

versus  k  corresponding  to  a  perfectly  conducting  waveguide  and 

2  2  2  2  2  1/2 

w  =  (k  c  +  ^  c  /Rq)  versus  k  corresponding  to  a  dielectric  loaded 
waveguide.  The  straight  line  w  *=  kB^c  represents  the  free-streaming 


mode.  In  a  range  of  physical  parameters,  the  mode 


kB,c  intersects 
0 


2  2  2  2  2  1/2 

uj  *  (k  c  +  £;  c  /Rq)  at  k  «  kp,  indicating  a  possible  mode  coupling. 
In  fact,  for  k  >  k^  in  Fig.  12,  the  phase  velocity  of  the  vacuum 
dielectric  mode  is  less  than  the  beam  velocity.  In  this  regard, 
we  expect  a  strong  Cherenkov  radiation'*'0 ’  near  the  intersection 

point  of  these  two  modes. 

In  order  to  investigate  stability  properties  of  the  space  charge 
wave  in  a  dielectric  loaded  waveguide,  it  is  necessary  to  numerically 
solve  Eq.  (18)  with  no  a  priori  assumption  that  the  beam  is  very  tenuous. 
However,  use  is  made  of  the  fact  that  the  Doppler-shifted  eigenf requency 


u  -  kB^c  is  well  removed  from  the  free-streaming  mode,  i.e.,  |u  -  kB^c|  << 
kB^c.  Evaluating  the  parameter  q  in  Eq.  (47)  and  the  wave  admittance 
F  in  Eq.  (19)  at  w  *  ojq  =  kB^c,  the  dispersion  relation  in  Eq.  (18) 
can  be  approximated  by 


(£0  +  6C  )  J1(C0  +  «C) 


where 


Fn  +  (DF/3w)  (w  -  kfi.c)  , 
0  w0  b 


F(w0,k)  , 


the  parameter  £q  *  C  (u)q)  is  defined  by 


w '  F°  ’ 


(56) 

(57) 

(58) 


and 
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6^  -  £  -  S 


0  • 


(59) 


Taylor  expanding  the  left-hand  side  of  Eq.  (56)  about  C  *  Cq  and 
making  use  of  Eq.  (58),  we  approximate  Eq.  (56)  by 


2  2 
5  “  *0 


1  + 


*o  +  Fo 


■r  OF/3w)  (u)  -  kB.c) 
2  w0  b 


(60) 


where  use  has  been  made  of  the  assumption  that  the  term  proportional 

to  (3F/3w)  in  Eq.  (60)  is  much  less  than  unity.  Substituting 
u0 

Eqs.  (14)  and  (20)  into  Eq .  (60),  we  finally  have  the  dispersion  relation, 


(j^(ui  -  kBbc  +  ikA/y^m)^ 


(‘ 


‘2-7'o 


(61) 


1  +  —5 - r  OF/Su })  (cu  -  k6,c) 


*0  +  Fo 

for  the  space  charge  wave  in  a  dielectric  loaded  waveguide. 

Defining  the  normalized  Doppler-shi f ted  eigenf requency  12  by 


12 


w  -  k6,  c 
_ b 

“pb/Yb 


(62) 


the  dispersion  relation  in  Eq.  (  61)  is  numerically  investigated  for 
a  broad  range  of  physical  parameters.  For  present  purposes,  to 
illustrate  the  mode  coupling  of  the  space  charge  wave  (o>  *  kB^c)  with 
the  vacuum  dielectric  mode,  shown  in  Fig.  13  are  plots  of  (a)  Fq  (solid 


curve)  and  (dashed  curve),  (b)  the  normalized  growth  rate 


lm!2 


and  (c)  Doppler-shifted  real  oscillation  frequency  ft =  Reft  versus  kR^ 
obtained  from  Eqs.  (57  ),  (58).  and  (61),  for  -y,  =  2,  e  ■  8,  R  /R  «=  0.8, 
Rw/Rc  =  0.8,  and  v  =  0.0025.  The  real  oscillation  frequency  in  Fig. 

13(c)  is  obtained  for  zero  axial  momentum  spread  (A  *  0).  In  Fig.  13(c), 
the  solid  curve  represents  the  unstable  mode  and  the  dashed  curves 
correspond  to  stable  oscillations.  Several  points  are  noteworthy  in 
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Fig.  13.  First,  the  maximum  growth  rate  for  instability  occurs  at 

2  1/2 

kR^  =  Yb(-t0)  corresponding  to  the  mode  coupling  point  in  Fig.  12. 

This  is  consistent  with  the  inductive  impedance  in  Sec.  III.  For 

example,  in  Fig.  13,  the  maximum  coupling  occurs  at  kR^  -  2  and 
2 

£q  “  -1.  Second,  the  maximum  growth  rate  is  order  of  the  beam  plasma 
frequency,  indicating  a  strong  instability.  In  this  regard,  this 
instability  can  be  utilized  to  generate  high  power  microwave.  Third, 
wavelength  of  the  microwave  radiation  generated  by  this  instability  can 
be  less  than  a  centimeter  for  a  subcentimeter  beam  radius.  Fourth, 
from  Fig.  13(c),  we  note  that  the  Doppler-shifted  real  frequency 
for  instability  is  negative,  thereby  implying  that  the  phase  velocity 
of  unstable  mode  is  less  than  the  beam  velocity.  We  therefore  conclude 
that  the  instability  mechanism  is  a  typical  Cherenkov  radiation. 

Finally,  the  growth  rate  and  bandwidth  of  instability  decrease  with 
increasing  value  of  the  axial  momentum  spread. 

The  dependence  of  stability  properties  on  the  ratio  Ro/R^  is 
further  illustrated  in  Fig.  14,  where  the  normalized  growth  rate  *  lmS2 
is  plotted  versus  kR^  for  A  =  0,  several  values  of  R^/R^,  and  parameters 
otherwise  identical  to  Fig.  13.  Obviously  from  Fig.  14,  we  note  that 
the  growth  rate  and  bandwidth  of  instability  increase  rapidly  as  the 
surface  of  dielectric  material  approaches  to  the  beam  surface  (Rn/R  ■*  1) 
for  a  given  beam  radius.  Shown  in  Fig.  15  is  plots  of  the  normalized 
growth  rate  versus  kR^  obtained  from  Eq .  (61)  for  y^  =  1.1547,  e  *  25, 
and  parameters  otherwise  identical  to  Fig.  13(b).  Even  for  a  moderate 
beam  energy,  the  growth  rate  of  instability  is  also  order  of  the  beam 
plasma  frequency.  However,  the  maximum  growth  rate  of  instability  is 
rapidly  decreasing  with  increasing  value  of  the  axial  momentum  spread. 
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Finally,  we  conclude  this  section  by  pointing  out  that  from  Figs. 

13(a)  and  (b),  the  maximum  coupling  of  instability  occurs  when  the  wave 
admittance  Fq  is  negative  [the  phase  angle  $  ■  ti  in  £q.  (23)]. 

In  this  regard,  the  root  £  can  be  a  pure  imaginary  value  (the  £r  «  0 
vertical  line  in  Fig.  1).  We  therefore  emphasize  that  the  wave 
admittance  F  in  a  dielectric  loaded  waveguide  can  be  perfectly 
Inductive  in  a  range  of  the  axial  wavenumber  k  corresponding  to  instability 
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VI.  CONCLUSIONS 


In  this  paper  we  have  investigated  properties  of  the  space  charge 

wave  in  a  solid  relativistic  electron  beam  propagating  in  a  cylindrical 

waveguide  with  an  arbitrary  impedance  Z.  The  perturbation  analysis 

was  carried  out  within  the  framework  of  the  linearized  Vlasov-Haxwell 

equations,  assuming  that  v/y^  <<  1.  The  formal  dispersion  relation  of 

the  space  charge  wave  for  azimuthally  symmetric  electromagnetic 

perturbations  (3/30  ■  0)  was  carried  out  in  Sec.  II,  including  the  important 

influence  of  finite  wall  impedance  Z.  In  Sec.  Ill,  the  resistive  wall 

instability  with  an  arbitrary  impedance  was  investigated,  showing  that 

an  inductive  impedance  wall  is  most  unstable.  Particularly,  the  maximum 

growth  rate  of  the  general  resistive  wall  instability  occurs  at  the 

2  2  2  2  2 

axial  wavenumber  k  satisfying  2k  RQ  =  Yb^i  ~  ^r^'  Properties  of  the 
space  charge  wave  in  a  perfectly  conducting  waveguide  was  investigated 
in  Sec.  IV.  It  was  found  that  the  space  charge  wave  admittance 
in  a  perfactly  conducting  waveguide  is  purely  capacitive,  thereby 
indicating  a  stable  propagation  of  the  electromagnetic  wave.  Moreover, 
we  obtained  the  condition  for  the  zero  phase  velocity  (w/k  *0),  in 
connection  with  collective  ion  acceleration.  Furthermore,  it  was  also 
shown  in  Sec.  IV  that  for  short  axial  wavelength  perturbations  (kR^/y^  >  10), 
the  eigenfunction  can  be  represented  by  a  Bessel  function.  Space 
charge  wave  properties  for  a  dielectric  loaded  waveguide  were  investigated 
in  Sec.  V.  It  was  found  that  a  strong  mode  coupling  between  the 
vacuum  dielectric  waveguide  and  beam  streaming  modes  occurs  in  a  range 
of  physical  parameters,  exhibiting  possibilities  of  a  strong  Cherenkov 
radiation.  The  maximum  growth  rate  of  instability  is  order  of  the 
beam  plasma  frequency.  In  this  regard,  the  Cherenkov  radiation  can 
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also  be  an  effective  means  to  produce  intense  high  power  microwave. 

The  wavelength  of  the  microwave  radiation  can  be  less  than  a  centimeter. 


ACKNOWLEDGMENTS 


This  research  was  supported  by  the  Independent  Research  Fund  at  the 


Naval  Surface  Weapons  Center. 


FIGURE  1  CONTOURS  OF  CONSTANT  PHASE  ANGLE  0  AND  MODULUS  f  [EQ.  (18)] 
IN  THE  COMPLEX  PLANE  $  =  (  £  r,  £  j)  FOR  n  =  1  RADIAL  MODE  NUMBER. 
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2 
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FIGURE  2  PLOT  OF  THE  NORMALIZED  GROWTH  RATE  Jl  j  VERSUS  y  [EQ.  (27)] 

FOR  £  2  >  £  2  AND  SEVERAL  VALUES  OF  5. 
i  r 
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8 


FIGURE  3  PLOT  OF  THE  PARAMETER  y  1  VERSUS  6  (SOLID  CURVE)  OBTAINED 

FROM  EQ.  (27)  FOR  £  2  >  £  2.  THE  DASHED  STRAIGHT  LINE  ISy'  = 
i  r  m 
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FIGURE  4  PLOT  OF  THE  NORMALIZED  GROWTH  RATE  ftj  VERSUS  y  [EQ.  (27)] 

FOR  £2  >£  2  AND  SEVERAL  VALUES  OF  6. 
r  i 
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8 


8 


4 
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FIGURE  5  PLOT  OF  THE  NORMALIZED  MAXIMUM  GROWTH  RATE  O .  AND 

CORRESPONDING  yc  VERSUS  6  [EQ.  (27)]  FOR  £2  >  t  2 
m  r  ’  i 
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FIGURE  6  PLOT  OF  THE  ADMITTANCE  F  =  Fr  VERSUS  £=  £  r  (DASHED  CURVES) 
OBTAINED  FROM  EQ.  (18)  AND  Fr  VERSUS  q  (SOLID  CURVES)  OBTAINED 
FROM  EQ.  (38)  FOR  R0/Rw  =  0.1,  0.5,  AND  0.9.  THE  HORIZONTAL  SCALE 
REPRESENTS  BOTH  £  r  AND  q. 
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0  1  2 


r/Ro 


FIGURE  8a  PLOT  OF  THE  EIGENFUNCTION  Ez  (r)  VERSUS  r/RQ  [EQ.  (17)]  FOR 
n  =  1,  R0/Rw  =  0.5,  AND  SEVERAL  VALUES  OF  q. 
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FIGURE  11a  PLOT  OF  THE  VACUUM  TM  MODE  DISPERSION  RELATION  IN  THE  PARAMETER 
SPACE  (w,k)  OBTAINED  FROM  EQS.  (48),  (50),  AND  (51)  FOR  Rw/Rc  =  0.8  AND 
SEVERAL  VALUES  OF  THE  DIELECTRIC  CONSTANT  e . 
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FIGURE  11b  PLOT  OF  THE  VACUUM  TM  MODE  DISPERSION  RELATION  IN  THE 

PARAMETER  SPACE  ( u>,k)  OBTAINED  FROM  EQS.  (48),  (50),  AND  (51) 
FOR  £  =  4  AND  SEVERAL  VALUES  OF  THE  RATIO  Rw/Rc. 
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FIGURE  12  SKETCH  OF  oj  =  (k2c2  +  /3 ^/R2)^  VERSUS  k  (CORRESPONDING  TO 
PERFECTLY  CONDUCTING  WAVEGUIDE)  AND  w  =  (k2c2  +  f  2c2/R2)Vi 
VERSUS  k  (CORRESPONDING  TO  AN  ARBITRARY  WALL  IMPEDANCE 
Z).  THE  STRAIGHT  LINE  u  =  k  0bc  IS  THE  FREE-STREAMING  MODE. 
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FIGURE  13a  PLOTS  OF  Fg  (SOLID  CURVE)  AND  {  J  (qaSHED  CURVE),  [EQS.  (57),  (58),  AND 
(61)]  FOR  7b  =  2,  e  =  8,  RQ/RW  =  0.8,  F?w/Rc  =  0.8,  AND  v  =  0.0025.  THE  REAL 
OSCILLATION  FREQUENCY  IS  OBTAINED  FOR  A  =  0. 


VERSUS  kR0  [EQS.  (57),  (58),  AND  (61 )]  FOR  7  b  =  2-  c  =  8-  R0/Rw  =  °-8- 


Rw/Rc  *  0.8,  AND  v  =  0.0025.  THE  REAL  OSCILLATION  FREQUENCY  IS 


OBTAINED  FOR  A  =  0. 
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FIGURE  14  PLOT  OF  THE  NORMALIZED  GROWTH  RATE  1m«  VERSUS  kRQ  [EQ.  (61)] 
FOR  A  =  0,  SEVERAL  VALUES  OF  RQ/RW,  AND  PARAMETERS  OTHERWISE 
IDENTICAL  TO  FIGURE  13. 
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